Comment on "Electrical transport in junctions between unconventional 
superconductors: Application of the Green function formalism" 
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Effect of surface pair breaking, entirely neglected 
by M. Samanta and S. Datta [ Phys.Rev. B 57, 10972 
(1998) ], is quite important in considering surface ( or 
interface ) quasiparticle bound states and associated 
characteristics of junctions involving unconventional 
superconductors. The whole class of bound states 
with nonzero energy is simply omitted within the 
framework of the approach, using uniform spatial pro- 
file of the order parameter up to the interface. The 
contribution of these bound states (as well as midgap 
states) to current-voltage characteristics of the SIS' 
tunnel junctions were studied for the first time in our 
earlier article. Dependence of midgap state contribu- 
tion to the Josephson critical current upon crystal 
to interface orientations is shown as well to be fairly 
sensitive to the effect of surface pair breaking. 

PACS numbers: 74.80.F, 74.50. +r 



In a recent article M. Samanta and S. Datta con- 
sidered theoretically electrical transport of junctions in- 
volving unconventional superconductors. In particular, 
they discussed contributions to junction properties from 
midgap surface states, that arise in d-wave superconduc- 
tors due to sign change of the order parameter. They 
pointed out that the effect of midgap states is most 
prominent for weakly coupled junctions (tunneling limit), 
concentrating in this respect mainly on the first order 
theory in transmission coefficient. In evaluating electric 
current across the junction the authors [Q neglected from 
the very beginning all the effects of surface pair breaking, 
considering order parameters from both sides of the junc- 
tion to be equal to their bulk values up to the junction 
barrier plane. In this Comment we demonstrate that the 
approximations jl) lead to incorrect results for uncon- 
ventional superconductors, since effects of surface pair 
breaking is of crucial importance for the I-V characteris- 
tics of tunnel junctions, especially due to the appearence 
of surface quasiparticle states with nonzero energy. Also 
we point out that correct theory for current-voltage char- 
acteristics of tunnel junctions involving anisotropically 
paired superconductors was developed for the first time 
in§. 

In contrast to s-wave isotropic superconductors, d- 
wave superconductors are known to be quite sensitive to 
any inhomogeneities (impurities, surfaces, interfaces). In 
particular, for many of crystal to surface orientations the 



order parameter turns out to be essentially suppressed on 
the tunnel barrier plane. Several important experimen- 
tal methods used for studying the anisotropic structure 
of the order parameter, for example, tunneling measure- 
ments are, in turn, fairly sensitive to the superconduct- 
ing properties close to the surface of the sample. The 
effects of anisotropic pairing on the tunneling density of 
states (the local quasiparticle spectrum at the surface), 
the Josephson and quasiparticle current of SIS' and SIN 
tunnel junctions were theoretically studied by taking ac- 
count of surface pair breaking and quasiparticle surface 
bound states in 

Our main assertion is that the whole class of surface 
quasiparticle states is omitted in JjJ due to the disregard 
the surface pair breaking there. Spatial profile of the or- 
der parameter suppressed near the surface, can be consid- 
ered as effective potential well for quasiparticles, which 
suffer Andreev reflection towards the surface in this case. 
Andreev reflection processes along with the conventional 
reflection from the surface, can result in forming quasi- 
particle bound states (Andreev bound states) localized 
near the surface within the characteristic length roughly 
of order of the superconducting coherence length. Quasi- 
particle surface bound states with nonzero energy are 
present in the case of surface pair breaking and do not 
exist for a uniform spatial profile of the order parameter. 
Only mid-gap surface states, having supersymmetric ori- 
gin, still exist for that uniform model profile. Thus, for 
the order parameter, which is independent of the spatial 
coordinate up to the surface or interface, one can find in 
the tunneling limit only peak at zero energy in the local 
density of states, while all nonzero peaks taking place in 
the presence of surface pair breaking turn out to disap- 
pear in the model. 

Bound states with nonzero energy result in the anoma- 
lies of current-voltage characteristics in the presence of 
externally applied voltage, described in and entirely 
omitted in Q). Positions and characteristics of those 
anomalies turn out to be associated with positions and 
types of extremal points of momentum dependence (dis- 
persion) of bound state energies. Possible characteristic 
points on the I-V curves for SIS' and SIN tunnel junctions 
are described in by using Eilenberger's equations 

for the quasiclassical Green's functions combined with 
corresponding boundary conditions and the microscopic 
expression for the tunnel current. 

Midgap states are dispersionless bound states, so that 
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their contributions to junction characteristics can dif- 
fer from the ones with nonzero energy. In particular, 
if midgap states take place only from one side of a 
SIS' tunnel junction, they contribute to the specific fea- 
tures of quasiparticle current (and not to ac Josephson 
effect) at the voltage values, determined by anomalies 
of the local density of states on the other bank of the 
junction. For example, in the case of an isotropic s- 
wave superconductor on the other bank, with an order 
parameter A, quasiparticle current has peaks at volt- 
age values eV = ±A. The peaks have inverse square- 
root "divergent" behavior in the vicinity of \eV\ = |A|: 
I w oc Q(\eV\ - \A\)/^\eV\ 2 - |AJf . This particular re- 
sult, obtained for the first time in g, was rederived in Jl]] 
in disregarding the surface pair breaking. There are var- 
ious reasons for the modification of this inverse square 
root singularity, at least sufficiently close to the point 
\eV\ = |A|, in particular, due to the surface roughness 
and bulk impurities. We do not discuss in the Com- 
ment this kind of effects, determining the hight of the 
peaks, since it is not discussed at all in Qj. We note that 
even in the absence of any surface pair breaking there are 
some additional peaks and jumps of the conductance of 
SIS' junction, involving anisotropically paired supercon- 
ductors. Positions of the specific features of the conduc- 
tance are determined by the extremal points of the sum 
of order parameters from both sides of the junction (and 
for the difference as well, although not for sufficiently 
low temperatures), taken for incoming and transmitted 
quasiparticle momenta M. These specific features are 
omitted in jlj as well. 

As it is known, for SIN tunnel junction of high qual- 
ity involving a d-wave superconductor, the zero-bias 
anomaly of the conductance takes place due to the 
midgap states. Let the normal metal suffer the super- 
conducting transition at lower temperature into a s-wave 
isotropic superconducting state. Then, according to the 
result of Ref. || just stated above, the zero-bias con- 
ductance peak should split into two peaks lying at the 
applied voltage eV = ±A(T). This splitting increases 
along with A(T) with further decrease of temperature. 
The behavior of I-V curves of this type for the tunnel 
junction between BiiSriC 'aC^Oa (hight-T c supercon- 
ductor) and Pb (s-wave superconductor) was recently 
observed experimentally in 0). There are several ex- 
perimental results, which could be interpreted as due to 
zero energy bound states [pl-fbif. At the same time no 
anomaly in the Josephson critical current due to midgap 
states is observed to our knowledge. Manifestations of 
bound states with nonzero energies are not noticed ex- 
perimentally up to now as well. Possibly, this is due to 
reasons resulting in brodening of Andreev bound states 
( e.g. due to interface roughness HQ ). 

In the presence of midgap states from both sides of 
tunnel junction, they cause the specific changes both in 
the quasiparticle current and in ac Josephson effect 



The effect of finite transmission of the barrier plane be- 
yond the tunneling limit may result in the shift of midgap 
states to nonzero values of energy ( though for particular 
case of "mirror" tunnel junctions and no phase difference 
beteen superconductors there is no shift for the midgap 
states due to this reason ). These "former midgap states" 
take place for the uniform model as well, in contrast to 
other interface bound states. 

Disregarding the surface pair breaking is the common 
feature of many articles, which consider surface (inter- 
face) bound states in <i-wave superconductors both in 
studying current-voltage characteristics |I],|l5|-|2"o|] and dc 
Josephson effect [^1-24 . It is worth noting, that the ef- 
fect of surface pair breaking can be of importance not 
only for studying the current-voltage characteristics but 
in considering the Josephson critical current as well. In 
particular, in the presence of essential suppression of the 
order parameter at the surface there are deviations of 
temperature dependence of the Josephson critical current 
I c from conventional Ambegaokar-Baratoff behavior both 
near T c |3| and at low temperatures || . For example, in 
the vicinity of T c the Josephson critical current turns out 
to be proportional to (T c — T) 2 for sufficiently small sur- 
face values of the order parameters on both banks of the 
junction. It is proportional to (T c — T) 3 / 2 if the order 
parameter vanishes only from one side of the tunnel bar- 
rier plane. Conventional Ambegaokar-Baratoff behavior 
I c oc (T c — T) holds near T c in the absence of any essential 
suppression of the order parameters at the interface. 

Since the low-temperature anomaly of the Josephson 
critical current is associated with the effect of midgap 
states Q, the influence of surface pair breaking on the 
characteristics of midgap states should be discussed in 
this context. The occurrence of the zero-energy peak in 
the tunneling density of states is unaffected by the self- 
consistency of the order parameter. However, disregard- 
ing surface pair breaking can result in essential overesti- 
mating the weight of the peak and, as a consequence, the 
Josephson critical current. Moreover, since surface pair- 
breaking is sensitive to the crystal to surface orientation, 
disregarding its effect results in qualitative changes in de- 
pendences of the peak height and I c upon the misorien- 
tation angles of superconductors from both sides of the 
junction. To illustrate the above assertion we consider 
the midgap state contribution to the Josephson critical 
current between two identical d-wave superconductors, as 
a function of the crystal orientation of superconducting 
electrodes with respect to the specularly reflecting tun- 
nel barrier plane. For the sake of simplicity we consider 
symmetric tunnel junction, when crystal axes have the 
same orientation in both electrodes. 



2 



3.5 



3.0 - 



2.5 - 



2.0 - 



1.5 - 



: . o 



0.5 - 



0.0 




5 10 15 20 25 30 35 40 45 

9 (degree) 

Fig. 1. Midgap state contribution to the Josephson critical 
current 1(9) as a function of misorientation angle 9 at temper- 
ature T = 0.1T C . Transmission coefficient D is taken oc cos 2 ip, 
where angle <p assigns the incoming momentum direction with 
respect to the interface normal. Critical current I CO nst(9) for 
the case of constant order parameter (dashed line) is com- 
pared with calculations using a self-consistently determined 
order parameter (lower solid line). Functions I aB if(ff) and 
Iconat(9) are normalized to the value I 3e if(9 = 45°). 

We use a cylindrical Fermi surface with an isotropic 
Fermi velocity and a constant normal state density of 
states. The calculation of the spatially dependent or- 
der parameter involves the pairing interaction, which we 
write as V{p f ,p' f ) = 2Vcos(20 - 29) cos(2</>' - 26). Here 
is the misorientation angle, that is the angle made by 
the interface normal with the principal axis of the <i-wave 
order parameter on each side of the junction. In Fig. 1 
we present the numerical results for the midgap states 
contribution to the Josephson critical current ( which 
dominates for sufficiently small temperatures ) , as a func- 
tion of the misorientation angle. The lower solid line de- 
scribes self-consistent results, while the dashed line cor- 
responds to neglecting the effect of pair breaking. The 
quantities are normalized to the self-consistent value of 
the Josephson critical current for 6 = 45°. The figure 
shows strong deviation of self-consistent angular depen- 
dence from the non-selfconsistent one. This reflects the 



fact that suppression of the order parameter near the 
tunnel barrier is essential for a wide range of crystalline 
orientations and sensitive to change of the orientation. 
The most significant deviation takes place for 9 = 45°, 
when the order parameter is completely suppressed near 
the boundary. For this particular orientation disregard- 
ing the surface pair breaking results in overestimation of 
midgap states contribution to the Josephson critical cur- 
rent more than in three times. Similar difference can be 
seen in Fig. 2 in Ref. ^5|. Our results demonstrate, in 
particular, the failure of the simple orientation depen- 
dence I c (midgap) oc sin(20i) sin(2# 2 ), obtained in [p3"P] 
on the basis of a uniform model order parameter and 
some additional approximations. One can show, how- 
ever, that signs of I c (midgap) and sin(20i) sin(2#2) coin- 
cide for the particular pairing potential considered. 

The authors of Ref. pi make also an accent on the 
fact that they do not use quasiclassical approximation in 
their theory. Extension of the theory beyond the quasi- 
classical approximation would be interesting both in the 
case of any new qualitative features of the phenomena in 
question, or due to noticeable quantitative corrections, 
containing powers of the parameter a/£, where a is the 
atomic scale. However, no new qualitative features of 
nonquasiclassical nature are found in [EJ. We believe, 
that taking account of the effect of surface pair break- 
ing within the quasiclassical theory, is much more impor- 
tant from the qualitative point of view, while quantita- 
tively disregarding this effect makes the extension of the 
approach beyond the quasiclassical approximation to be 
senseless. 
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